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AERONAUTIC  SYMBOLS 

I.  FUNDAMENTAL  AND  DERIVED  UNITS 


Metric 

English  j 

Symbol 

Unit 

Abbrevia¬ 

tion 

Unit 

Abbrevia¬ 

tion 

Length _ 

l 

meter _ _  _  __ 

m 

foot  (or  mile)  _ _ 

ft.  (or  mi.) 

Time _ 

t 

second..  ______ 

s 

second  (or  hour)  _ 

sec.  (or  hr.) 

Force _ 

F 

weight  of  1  kilogram. 

kg 

weight  of  1  pound  _  __ 

lb. 

Power _ 

P 

horsepower  (metric)  __ 

horsepower _ 

hp. 

Speed _  _ 

V 

/kilometers  per  hour 

k.p.h. 

miles  per  hour _ 

m.p.h. 

tmeters  per  second _ 

m.p.s. 

i 

feet  per  second _ 

f.p.s. 

2.  GENERAL  SYMBOLS 


Weight  =  mg 

Standard  acceleration  of  gravity  =  9.80665 
m/s2  or  32.1740  ft./sec.2 

Mass  =  T 

3 

Moment  of  inertia  =  mk2.  (Indicate  axis  of 
radius  of  gyration  k  by  proper  subscript.) 
Coefficient  of  viscosity 


v,  Kinematic  viscosity 

p,  Density  (mass  per  unit  volume) 

Standard  density  of  dry  air,  0.12497  kg-m“4-s2  at 
15°  C.  and  760  mm;  or  0.002378  ib.-ft.“4  sec.2 
Specific  weight  of  “standard”  air,  1.2255  kg/m3  or 
0.07651  Ib./cuit. 


3.  AERODYNAMIC  SYMBOLS 


Area  of  wing 
Gap 
Span 
Chord 

Aspect  ratio 

True  air  speed. 

Dynamic  pressure  =^pV2 

Lift,  absolute  coefficient  CL  =  ~ 

qS 

Drag,  absolute  coefficient  CD  =  Q 

q& 

Profile  drag,  absolute  coefficient  CDo  =  ^ 
Induced  drag,  absolute  coefficient 

1  q& 

Parasite  drag,  absolute  coefficient  Cn 

p  q$ 

Cross-wind  force,  absolute  coefficient 

gS 

Resultant  force 


Angle  of  setting  of  wings  (relative  to  thrust 
line) 

Angle  of  stabilizer  setting  (relative  to  thrust 
line) 

Resultant  moment 
Resultant  angular  velocity 

Reynolds  Number,  where  l  is  a  linear  dimension 
(e.g.,  for  a  model  airfoil  3  in.  chord,  100 
m.p.h.  normal  pressure  at  15°  C.,  the  cor¬ 
responding  number  is  234,000;  or  for  a  model 
of  10  cm  chord,  40  m.p.s.  the  corresponding 
number  is  274,000) 

Center-of-pressure  coefficient  (ratio  of  distance 
of  c.'p.  from  leading  edge  to  chord  length) 
Angle  of  attack 
Angle  of  downwash 
Angle  of  attack,  infinite  aspect  ratio 
Angle  of  attack,  induced 

Angle  of  attack,  absolute  (measured  from  zero- 
lift  position) 

Flight-path  angle 
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A  SIMPLIFIED  APPLICATION  OF  THE  METHOD  OF  OPERATORS  TO  THE 
CALCULATION  OF  DISTURBED  MOTIONS  OF  AN  AIRPLANE 

By  Robert  T.  Jones 


SUMMARY 

A  simplified  treatment  of  the  application  of  Heavisides 
operational  methods  to  problems  of  airplane  dynamics  is 
given.  Certain  graphical  methods  and  logarithmic  for¬ 
mulas  that  lessen  the  amount  of  computation  involved  are 
explained . 

The  problem  of  representing  a  gust  disturbance  or 
control  manipulation  is  taken  up  and  it  is  pointed  out 
that  in  certain  cases  arbitrary  control  manipulations  may 
be  dealt  with  as  though  they  imposed  specific  constraints 
on  the  airplane ,  thus  avoiding  the  necessity  of  any  inte¬ 
gration  whatever. 

The  application  of  the  calculations  described  in  the  text 
is  illustrated  by  several  examples  chosen  to  show  the  use 
of  the  methods  and  the  practicability  of  the  graphical  and 
logarithmic  comp utations  described . 

INTRODUCTION 

The  theory  of  airplane  dynamics  in  its  present  form 
is  due  mainly  to  the  original  researches  of  Lanchester 
and  Bryan  on  the  stability  of  airplanes.  Later  inves¬ 
tigators,  notably  Bairstow  and  Wilson  (reference  1), 
applied  and  extended  the  original  conceptions  of  the 
theory.  Bryant  and  Williams  (reference  2)  have  re¬ 
cently  shown  how  the  operational  mathematics  of 
Heaviside  may  be  used  in  applying  the  theory  to  prob¬ 
lems  of  the  disturbed  motions  of  airplanes. 

Although  the  calculation  of  disturbed  motions  of 
aircraft  is  important  in  problems  of  flight  safety,  little 
experience  has  been  gained  in  the  practical  application 
of  the  theory  owing  to  its  mathematical  complexity. 
The  present  paper  gives  the  results  of  researches  in  the 
mathematical  application  of  the  theory.  It  has  been 
found,  as  suggested  by  Bryant  and  Williams,  that  the 
Heaviside  method  affords  the  simplest  and  most  direct 
solution  of  these  problems.  In  order  to  bring  out  the 
advantages  of  this  method,  a  treatment  of  its  applica¬ 
tion  is  given  and  certain  formulas  and  graphical  con¬ 
structions  are  explained  that  make  the  calculations 
easier. 

In  their  usual  form,  problems  of  airplane  dynamics 
depend  for  solution  on  -the  integration  of  simultaneous 
linear  differential  equations.  Methods  for  the  integra¬ 


tion  of  such  equations  are  given  by  Wilson  and  Bouth 
(references  1  and  3)  and  in  mathematical  textbooks. 
The  problems  met  in  airplane  dynamics  are  often  more 
complex  than  the  examples  treated  in  textbooks  and, 
when  an  attempt  is  made  to  apply  the  given  methods 
to  their  solution,  difficulties  of  computation  arise. 

In  view  of  the  importance  of  investigating  these 
problems  and  since  their  solutions  involve  lengthy  cal¬ 
culations,  it  is  desirable  that  as  many  mathematical 
simplifications  as  possible  be  employed.  Heaviside’s 
method  gives  such  a  simplification,  the  solution  of  the 
differential  equations  being  accomplished  symbolic  all}7 
by  a  single  “expansion  theorem.” 

THE  DIFFERENTIAL  EQUATIONS  FOR  THE  DISTURBED 
MOTIONS 

An  airplane  in  uniform  flight  may  be  thought  of  as  a 
free  rigid  body  in  equilibrium.  Deviations  of  the  air¬ 
plane  from  this  equilibrium  condition  may  be  caused 
by  reactions  due  to  control  movement,  gustiness  in  the 
air,  or  by  some  influence  such  as  the  stopping  of  an 
engine.  The  motions  of  the  airplane  following  such  a 
disturbance  may  be  calculated  if  the  momentary 
accelerations  or  forces  are  known.  It  is  obvious  that 
this  computation  may  be  performed  by  taking  small 
intervals  of  the  time  and  caculating  the  velocities  and 
displacements  generated  by  the  known  accelerations 
step  by  step,  assuming  the  accelerations  momentarily 
constant. 

The  component  linear  and  angular  motions  of  the 
airplane  in  its  deviations  from  equilibrium  are  given 
exact  definition  by  constructing  a  set  of  axes  rigidly 
fixed  in  the  machine  and  considering  its  motions  as 
being  those  of  the  axes  themselves.  The  motions 
spoken  of  are  then  velocities  and  displacements  of  the 
airplane  axes  relative  to  the  earth  or  the  air.  When 
the  airplane  is  in  steady  flight,  it  maintains  a  certain 
equilibrium  attitude  with  respect  to  the  air  and  to 
the  earth.  Thus  for  climbing  flight  at  a  given  engine 
speed  a  definite  angle  of  attack  and  a  definite  angle  of 
pitch  must  be  preserved.  Deviations  from  equilibrium 
in  either  sense  will  introduce  reactions;  hence  motions 
of  the  airplane  axes  relative  to  both  air  and  earth 
must  be  considered. 
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The  aerodynamic  reactions  to  the  motions  arise 
from  changed  relative  air  velocities  over  the  different 
parts  of  the  airplane.  The  calculation  or  measure¬ 
ment  of  these  component  aerodynamic  reactions  leads 
to  quantities  known  as  “resistance  derivatives”  or 
“stability  derivatives,”  which  are  taken  as  constant 
factors  of  proportionality  between  the  reactions  and 
the  velocities  or  displacements  of  the  motions.  For  a 
more  detailed  exposition  of  the  concept  of  stability 
derivatives,  the  reader  is  referred  to  standard  text¬ 
books  on  aeronautics. 

On  account  of  the  bilateral  symmetry  of  the  airplane 
it  is  customary  to  divide  the  motions  into  two  inde¬ 
pendent  groups,  the  lateral  and  the  longitudinal,  each 
consisting  of  three  degrees  of  freedom: 

f  Rolling. 

(A )  La  ter  al  mo  tions _ Y  awing . 

[Sideslipping. 

Pitching. 

(B)  Longitudinal  motions _ Vertical  translation. 

Forward  translation. 
Presumably  the  reactions  to  small  increments  of 
longitudinal  speed  or  displacement  do  not  sensibly 
influence  the  lateral  motions  and  the  two  groups  may 
be  independently  treated.  In  order  to  illustrate  the 
calculation  of  the  history  of  a  motion  due  to  a  given 
disturbance,  examples  of  lateral  motions  are  chosen 
although  the  methods  used  are  equally  applicable  to 
any  set  of  degrees  of  freedom  of  the  airplane.  The 
quantities  that  arise  in  the  consideration  of  the  lateral 
motions  are  defined  in  the  following  table: 

Velocities  and  displacements  of  airplane  axes: 

U(t,  equilibrium  flight  velocity  along  Ar  axis. 

/*,  component  of  flight  velocity  along  F  axis 
(sideslipping). 

p,  component  of  angular  velocity  about  X  axis 
(rolling). 

ry  component  of  angular  velocity  about  Z  axis 
(yawing). 

<p}  angle  of  bank  (relative  to  gravity). 

Forces  and  moments  resolved  along  airplane  axes: 

F,  component  of  force  along  Y  axis. 

L ,  component  of  moment  about  X  axis  (rolling 
moment). 

N,  component  of  moment  about  7  axis  (yawing 
moment). 

Accelerations  of  airplane: 

Yo—Y/m  (force  per  unit  mass). 

L0=L/mkx 2  (moment  per  unit  moment  of  inertia). 

N0=N/mkz2  (moment  per  unit  moment  of  inertia). 

Gust  velocities  resolved  along  airplane  axes: 


Note. — The  signs  of  the  gust  velocities  are  so 
chosen  that  a  positive  gust  produces  the  same  aero¬ 
dynamic  reaction  on  the  airplane  as  a  positive  velocity 
of  the  airplane  in  still  air.  The  resolution  of  gust 
velocities  along  the  moving  axes  is  exact  only  to  the 
first  order  of  the  small  quantities  involved. 

Airplane  characteristics  used  as  parameters: 

Yv 

Stability  derivatives  in  terms  of  accelera¬ 
tor  tions  of  airplane,  thus: 

Lp  „  dYI 

do  / 

Lv 

L, 


With  the  definition  of  the  component  motions  that 
are  to  be  considered,  the  stability  derivatives  will  be 
of  the  form: 

dL  dN  dY  , 

^  ’  e  t  e . 

dp  do  d  r 

where  Z,  A7,  F,  respectively,  are  the  rolling  moment, 
the  yawing  moment,  and  the  sidewise  force,  as  they  are 
customarily  defined. 

It  has  been  found  convenient  to  transform  all 
stability  derivatives  and  disturbing  effects  into  terms 
of  accelerations  of  the  airplane  rather  than  retaining 
them  as  moments  and  forces.  This  transformation 
is  accomplished  by  dividing  out  the  appropriate 
moments  of  inertia  and  the  mass  of  the  machine. 

For  example,  jmkx2  may  be  written  simply  as  Lp) 

similarly  jmkz2=Nv  and  ^-j  m  =  Yr- 

If  the  flight  path  is  assumed  to  be  horizontal  (or 
nearly  so)  and  the  main  forward  velocity  £70  to  be 
substantially  constant,  the  equations  of  motion  in  a 
lateral  disturbance  may  be  written: 

(In  sideslipping)  =  g# — r  Un + v  Y„ + r Y,  +  Y0  j 


(In  rolling) 


- — vL  v + pLp + rLT +Z0 


(In  yawing)  -^t^vNv+pNp + rNf + N0 


- o }  component  of  gust  velocity  directed  along  F  In  these  equations  the  terms  F0,  L0,  and  Ar0  represent 
axis.  known  disturbing  or  controlling  accelerations,  assumed 

r0,  component  of  angular  velocity  of  gust  about  Z  to  be  given  as  functions  of  the  time  t.  In  the  first 
axis.  equation  the  terms  g<p  and  —rU0  are,  respectively,  the 

Pot  component  of  angular  velocity  of  gust  about  X  accelerations  due  to  gravity  and  to  the  rotation  of  the 
axis.  moving  axes.  Since  the  axes  chosen  will  ordinarily  lie 
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near  the  axes  of  the  principal  moments  of  inertia  of 
the  airplane,  terms  involving  the  products  of  inertia 
have  been  neglected. 

INTEGRATION  OF  EQUATIONS  FOR  VELOCITIES  AND 
DISPLACEMENTS 

As  previously  mentioned,  equations  (1)  may  be  inte¬ 
grated  by  taking  small  intervals  of  the  time  and  calcu¬ 
lating  the  velocities,  and  finally  the  displacements,  by 

dv 

assuming  the  accelerations  -1^1  etc.,  to  be  momen¬ 
tarily  constant.  Although  this  method  is  sometimes 
useful,  it  naturally  leads  to  extensive  numerical  work. 
The  operational  mathematics  of  Heaviside  appear  to 
offer  the  most  promising  means  of  performing  these 
integrations. 

The  first  step  in  integrating  the  equations  of  motion 
by  the  operational  method  is  to  replace  the  symbol 

4-  by  the  so-called  “differential  operator”  D ,  which  is 
dt 

to  be  treated  as  though  it  were  an  ordinary  algebraic 
quantity;  the  equations  are  then  rearranged  with  the 
known  disturbance  effects  on  the  right-hand  side: 

(D-Yv)v-g<p+  (U0-  Yr)r—  F0 

—Lvv+  (D—Lp)p—Lrr=L0  (2) 


In  the  calculation  of  any  of  the  velocity  components 
the  same  denominator  appears;  if  this  determinant  is 


denoted  by  F(D),  the  forms  of  these  components  are: 

„  _Jn(D) y  i  f  12(D)  t  .  j\ 3 CP)  at 

FID r  0+ F(D)U^F(D)  0 

(5) 

/21  (D)  y  I J 22(D)  j  .  /23(D)  w 
p~ F(D)lo+F(D)  F(D)  0 

etc. 

Thus  far  the  solution  of  the  equations  of  motion  has 
progressed  simply  on  algebraic  grounds,  the  required 
quantities  (v,  p,  etc.)  having  been  found  explicitly  in 
terms  of  the  symbol  D.  The  symbol  D  was  defined  as 
the  operation  of  derivation  with  respect  to  the  time  t , 
expressed  by  writing 

t  1  d 


The  terms  of  the  solution  j(D)jF(D)  indicate  that  the 
formal  operations  are  to  be  performed  on  whatever 
functions  follow  them  as  factors.  Since  they  contain 
the  symbol  D  in  their  denominators,  it  becomes  neces¬ 
sary  to  define  the  operation  indicated  by  1  ID  or  D~K 
As  D  is  an  operation  and  not  a  number,  its  reciprocal  is 
defined  as  the  inverse  of  the  operation  of  differentiation, 
rather  than  as  the  derivative  itself  divided  into  1. 
The  inverse  of  differentiation  is  integration:  thus, 


—Nvv—Npp+  (D—-NT)r—N^ 

Since  D<p=p ,  the  first  equation  may  be  operated  on 
throughout  by  Z>,  reducing  all  to  the  same  variables 
(v,  p ,  r): 

D(D-Y9)v-gp+D(UQ-Yr)r=DY0 

— L  Vv  T  (D  Lp)p  Lrr — L0  (2  a) 

-Nvv-Npp+(D-NT)r=N0 

With  the  equations  in  this  form,  they  may  be  solved 
for  v,  p,  or  r  by  ordinary  algebraic  means;  thus, 

DY0  -g  D  (  Uq  Y r) 

U  (D—Lp)  -Lr 

N0  —Np  (D—Nr) 

- t - (3) 

D(D—  Yv)  -g  D(U0-Yr) 

—Lv  (D—Lp)  - Lr 

~NV  -Np  (D-Nr) 

The  expansion  of  the  determinant  of  the  numerator 
in  terms  of  minors  results  in: 


ZF1-/  .  .  .  dt 

The  operations  indicated  by  .  the  ratios  of  polynomials 
in  D  that  occur  in  the  terms  of  our  solution  then  consist 
of  a  succession  of  differentiations,  [f(D)]f  and  a  succes¬ 
sion  of  integrations,  [F(D)}~1.  It  is  clear  that  the 
nature  of  the  problems  at  hand  requires  that  the  result¬ 
ant  of  these  operations  be  an  integration,  which  is 
shown  by  the  fact  that  the  polynomial  F(D)  is  invar¬ 
iably  of  higher  degree  in  D  than  any  of  the  polynomials 
f(D). 

THE  EXPANSION  EQUATION 

By  treating  the  disturbances  (such  as  Y0j  N0 ,  L0)  as 
discontinuous  functions  of  the  time,  Heaviside  obtained 
solutions  of  equations  similar  to  the  foregoing  by  a 
simple  theorem.  The  substitution  of  Y0  into  Heavi- 
side’s  theorem  results  in 


where  the  \Js  are  the  roots  of  the  polynomial  equation 
F(D)  =  0.  This  polynomial,  F(D)  —  0,  is  used  in  the 
study  of  the  stability  of  motion,  being  called  the 
“stability  equation.”  Its  roots,  X2  .  .  .  X*,  give  an 


|  (D-L,) 

-Lt  I 

I  D(U0-Yr ) 

-0 

-0 

D(U0-Yr ) 

\DYn+ 

Lq  + 

~Lr 

1  -A7, 

(D  —  Nr)  1 

1  (D-Nr) 
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indication  of  the  natural  tendencies  of  an  airplane’s 
motion  and  arc  used  in  the  definition  of  stability. 

In  order  to  apply  the  foregoing  theorem  to  the 
integration  of  equations  of  airplane  motion  it  is  neces¬ 
sary  to  assume  that  the  disturbance  terms  (Y0,  AT0, 
etc.)  due  to  the  control  or  gust  in  question  arc  instantly 
applied  at  the  assumed  origin  of  the  time  scale  (£~0) 
and  remain  constant  thereafter.  In  the  general  case 
the  disturbance  terms  in  the  equations  of  motion  cannot, 
be  thus  represented  as  remaining  constant  although  in 
practical  problems  they  may  almost  invariably  be 
represented  by  means  of  functions  of  the  form  Y0  ent. 
The  interpretation  of  Heaviside’s  theorem  equation  (6) 
when  this  form  of  function  is  used  is  (see  reference  2): 


When  dealing  with  variable  disturbance  terms,  it  is 
important  to  note  that  a  discontinuity  of  the  function 
representing  the  disturbance  at  t=0  is  implied  as  in 
equation  (6). 

By  the  substitution  of  (in)  for  n  in  equation  (7), 
expressions  that  can  be  used  when  the  disturbances 


Figure  l. — -Map  of  polynomial.  F(D)  =  Di-j-(iD>-\-bD2-{-cD-\-d  near  zero. 
F(D)->F(\)~0  when  D-tndzib. 

are  represented  by  forms  involving  sin  nt  or  cos  nt  are 
obtained 


= Y(jrn (cos  nt+i  sin  nt) +E(^)r%eX1](7a) 

X 

f  (vYl) 

If  —A+iB,  then  the  expressions  for  the  sine  and 
cosine  forms  separately  become 

(m^r  ^  . 

F(D)  ^  0  sm 

=  1  [a  Sin  nt+B  cos  nt + ^(xq^)/^(x)  ® 

,m)v 

F(L)')  ^  0  ^ 


F0  A  cos  nt~B  sin  nt- 


ZjU !+»>>'(x)  J 


These  latter  forms  are  particularly  useful  because 
almost  any  arbitrary  variation  of  gust  or  control  may 
by  expressed  as  a  sum  of  sine  or  cosine  terms.  Thus  if 

Y0~Ki  sin  nAFKn  sin  n2t-\-  .  .  .  +  etc. 


mv 

F(D) 1 0 


Ml 

F(D) 


Ki  sin  n-it 


m 

F(D) 


Ko  sin  n2t 


.  .  +  etc. 

(10) 


Each  of  these  terms  may  be  evaluated  by  equation  (8). 


SOLUTION  OF  OPERATIONAL  EQUATIONS 


FINDING  THE  ROOTS  OF  THE  EQUATION  F(Z))=0 

The  expansion  equations  given  for  the  forms 
f(D)/F(D)  require  the  roots  of  the  complementary 
equations  F(D)  =  0  for  their  solution.  In  cases  of 
airplane  motions  this  equation  is  normally  of  the  fourth 
degree  in  D;  hence  it  is  not  practicable  to  find  the  roots 
directly.  Although  a  number  of  methods  for  approxi¬ 
mating  the  roots  of  such  equations  have  been  devised, 
the  most  direct  way  is  to  draw  a  curve  of  the  function 
F(D)  against  D,  locating  the  real  roots  as  the  points 
crossing  the  I)  axis.  Usually  in  equations  of  this 
type  near  roots  may  be  isolated  by  separating  the 
equation  into  two  parts.  Thus,  if 


F(D)  =  IF + a  JY + bD1 +cD+d=Q  (11) 


there  will  usually  be  a  large  real  root  near  Z>4=  —aU) 

or  D——a,  and  a  small  one  near  D=-—-  This 

c 

division  follows  from  the  consideration  that  large  roots 
are  more  dependent  on  the  coefficients  of  the  higher 
powers  of  D  and  small  roots,  on  the  lower  powers. 

If  the  natural  motion  of  the  airplane  contains  oscilla¬ 
tory  components,  as  it  usually  does,  there  will  be  pairs 
of  conjugate  complex  roots  of  the  polynomial  F(D)  —  0 
in  addition  to  the  real  roots.  The  determination  of 
these  roots  is  naturally  more  difficult,  although  if  real 
roots  have  been  previously  found  they  may  be  used 
to  reduce  the  degree  of  the  equation  by  synthetic  divi¬ 
sion  and  the  determination  of  further  roots  will  be¬ 
come  progressively  easier.  Complex  roots  of  such  an 
equation  may  be  directly  found  by  plotting  a  map  of 
the  polynomial  F(D)  for  various  values  of  D  using  the 
coordinates  D=xJriy  and  finding  the  zero  point,  or 
root,  by  interpolation,  as  is  shown  in  figure  1.  If  a 
very  accurate  value  of  the  root  is  required  it  may  be 
convenient  to  plot  the  region  of  F(D)  near  the  origin 
to  a  magnified  scale.  Since  the  polynomial  is  what 
is  known  as  an  “analytic  function”  (reference  4), 


dF(D)  .dF(D) 
dx  1  dy 


(12) 


and  the  map  in  its  smallest  parts  will  consist  of  squares. 
In  this  way  a  more  accurate  interpolation  may  be 
made  or  a  process  analogous  to  Newton’s  method  may 
be  applied. 

It  will  be  found  most  convenient  to  calculate  the 
various  values  of  F(D)  by  means  of  a  vector  diagram 


the  calculation  of  disturbed  motions  of  an  airplane  o 

as  shown  in  figure  2.  If  trial  values  of  D  are  expressed  By  a  separation  of  the  two  components  of  the  imagi- 
in  the  form  R  (cos  6+i  sin  6)  or  Reie,  vectors  repre-  nary  root  \=a+ib,  the  logarithmic  formula  may  be 
senting  each  of  the  terms  of  the  polynomial  may  be  reduced  to 
simply  calculated.  The  problem  is  to  make  all  terms  l(a-\-jK) 

of  the  polynomial  balance  each  other  and  it  is  readily  l°g/a_i -ifyp' (a+ib)6^1^1^  (a~\~ib)tYlog  ^o+^o 

seen  how  this  may  be  accomplished  by  varying  6  to 

change  the  relative  inclinations  of  the  vectors  and  by  =  (log  R0+at)  +i(d0+bt)  (17) 

varying  R  to  change  their  relative  lengths.  The  Tbc  pna|  formu]a)  where  X, =a+ib,  then  becomes 
advantage  of  this  method  is  that  it  enables  a  close 

approximation  of  the  value  of  a  root  with  a  minimum  f(\ Q  g  (18) 

number  of  trials,  the  diagram  making  apparent  how  i) 

nearly  all  the  vectors  cancel  each  other.  OTj  by  De  Moivre’s  theorem, 

SOLUTION  OF  EXPANSION  EQUATIONS  „  . 

The  numerical  operations  indicated  in  the  expansion  \F' (xF  sm  (19) 

equations  (6)  to  (9)  call  for  calculations  with  complex  ,  ...  . 

numbers  (i.  e„  roots  of  F(Z»  =  0).  A  great  deal  of  the  The  points  thus  plotted  will  lie  on  a  logarithmic 


=R0ealel<-l’,^l,<>) 


nearly  all  the  vectors  cancel  each  other.  OTj  by  De  Moivre’s  theorem, 

SOLUTION  OF  EXPANSION  EQUATIONS  . 

The  numerical  operations  indicated  in  the  expansion  \F' (\fU  =K(>e"^C'°S  sm  (1!}) 

equations  (6)  to  (9)  call  for  calculations  with  complex  ,  ...  . 

numbers  (i.  e.,  roots  of  F(D)  =  0).  A  great  deal  of  the  The  points  thus  plotted  will  he  on  a  logarithmic 
labor  involved  in  these  computations  may  be  saved  by  spiral  (fig.  3);  the  deviation  of  this  spiral  from  a  circle 
the  use  of  graphical  and  logarithmic  methods.  CZ>,  i 

Thus,  if  it  is  desired  to  calculate  values  of  the  com-  '' 

plex  terms  occurring  in  equation  (6),  the  logarithmic  ^  \  y 

formula  cD+bDU''''  '  D* 


log^U^^M+log/fXd-logX.-logAAX,)  (13)  \  fV^\ 

\ _ ’ _ \ - + 

is  used.  For  the  purpose  of  calculating  these  loga-  \  \  3  4 

rithms,  it  is  convenient  to  express  the  complex  numbers  \  s'  3o)*\  /*'aD*+I>* 

(Xi,  /(Xi),  etc.)  as  vectors  of  radius  R  and  angle  9,  1  \.  / 

writing,  for  example,  1  .  aD3 

'K\  =  Cb ib =  R 1  (cos  #i  +  l  sin  0])  =R,e‘°i  (14)  FiauRE  2.— Graphical  method  of  locating  values  of  F(D)  neat  zero,  where 

by  De  Moivre’s  formula.  '  _ 

A  complex  term  of  equation  (6)  may  then  be  written,  gbows  tbe  jnfluence  of  damping  on  the  natural  motion 

/(X ,)  ^l(_  R2ea2  _x,  t  of  the  airplane. 

XiA'(Xi)e  R^eit-lR',fii,'i'  The  summation  indicated  in  equation  (6)  calls  for 

^  the  plotting  of  such  a  logarithmic  spiral  for  each  of  the 

= tt+i(b2—e1—d3)  (15)  complex  roots.  Since  these  roots  always  occur  in  con¬ 

jugate  pairs,  the  calculation  may  be  carried  out  for  one 
=A0ex,'+i9o  0f  sucb  a  pajr  an(j  a  spiral  calculated  for  the  second 

/(X  )  would  be  exactly  conjugate  to  the  first.  Thus,  it  is  only 

log  <eX|<— Xd+log  R0-\-id0  (16)  necessary  to  perform  the  foregoing  calculations  for  one 

1  1  root  of  each  pair,  the  summations  indicated  in  the 

and  the  resultant  logarithm  may  be  plotted  as  a  straight  e nuations  being  carried  out  in  effect  bv  merely  doubling 
line  Xif+ constant,  which  is  then  divided  or  extended  to  ^  abgcissas  o{  tbe  points  of  one  0f  the  conjugate 
represent  any  division  or  extension  of  the  time  t  over  -rajs  jf  ^  —a+ib  and  X2=a—  ib,  this  summation 
which  the  calculation  is  made.  (See  fig.  3.)  The  final  ,  ... 


aD* 


Figure  2. — Graphical  method  of  locating  values  of  F(D)  near  zero,  where 
F(D)=D*-\-aJ)3+bD?+cD±d.  Z>i=a:i+i3/i=I?ie^i=I3i(cos  9i+i  sin  0i) 


R2ei02 

~  Riei$iR^eie3 

^Rjtf'‘+' 

=R0e>"t+i80 


-eX|,+i(^— 8i— 83) 


«xji_  x^+log  Rq-^iSq 


vectors  will  represent  the  complex  values  of 

/(A) 

AF'(A)e 

and  it  is  seen  that  the  ordinates  of  the  points  of  the  line 
Xi<  +  constant  are  the  angles  of  these  final  vectors  while 
the  abscissas  are  the  logarithms  of  their  radii. 


A2 

^jxFAx)eXi=2^°ea<C0S  ^tJr6«)  (2°) 

>>1 

The  formulas  for  the  integration  of  terms  containing 
sin  nt  and  cos  nt  may  be  put  into  a  more  convenient 
form  for  the  graphical  or  logarithmic  calculations,  i.  e~, 


Y£[jjY0  sin  nt=  Y0  imaginary  coordinate  +  p,  (X)« 


Y0  cos  nt=  Y0  real  coordinate  of  p^J^eintYj 


-n2F(X)e 


6 


REPORT  NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


•III  these  forms  the  graphical  construction  of  the  terms 

T(m)€*nt  Proceecls  aI°ng  the  same  lines  as  that  of  the 

terms  involving  complex  roots  X.  Here  the  re¬ 
sulting  diagrams  will  be  circles,  divided  into  equal 
angles  as  nt  may  be  divided.  In  case  X  is  complex  the 
plot  of  an  ext  term  will  be  a  logarithmic  spiral  as  before 
and  it  is  important  to  remember  that  the  summation 


For  the  velocities  of  an  assumed  gust,  forms  involv¬ 
ing  enl  are  useful.  Thus,  if  the  gust  is  considered  a 
“transient”  one,  disappearing  rapidly  from  an  arbi¬ 
trary  initial  value,  the  form  (A)  (fig.  4) 

v0=Vie-nt  (24) 

may  be  used.  Here  i\  is  the  initial  value  and  -n  is 
chosen  to  make  the  gust  diminish  in  any  required  way. 


Figure  3.— Graphical  construction  of]Cw7Xc: 

X  1  (X) 


Log  /(X)-  log  X— log  F'{\)  =  log  i?n-h'0o 


over  eacn  or  a  pair  ot  conjugate  roots  is  accomplished 
by  doubling  the  abscissas  of  the  spiral  obtained  for  one. 

WAYS  OF  REPRESENTING  GUSTS  AND  CONTROL 
MANIPULATIONS 

GUST  DISTURBANCES 

If  the  disturbances  to  be  considered  are  due  to  gusts, 
the  terms  Y0,  L0,  etc.  of  equations  (1)  will  be  of  the  form 

I  (J  ~  2-'0  I  Y r 

(23) 

A  o — VqNv + V{\Nv + TqNt 

where  r0,  p0,  and  r0  are  the  component  velocities  of  the 
gust,  which  may  vary  with  the  time.  As  given,  the  ref¬ 
erence  system  for  specifying  these  gusts  has  been  chosen 
so  that  a  positive  gust  velocity  may  be  considered  as 
producing  the  same  aerodynamic  reaction  on  the  air¬ 
plane  as  a  positive  airplane  velocity  in  still  air.  All 
such  gusts  must  be  assumed  to  be  moderate  so  that 
second-order  effects  may  be  neglected.  (See  refer¬ 
ence  1.) 


If  the  gust  is  to  be  made  to  start  from  an  initial  value 
of  zero  and  to  persist  with  the  time,  the  form  (B) 

v0=vl(l—e-nt)  (25) 

may  be  used.  (See  reference  1.) 

For  the  purpose  of  representing  gusts  that  arise  with 
any  degree  of  sharpness  from  zero  velocity  to  a  given 
peak  value  and  then  diminish,  the  form  (C)  may  be 
used: 

v0=Ke~n((  (26) 

The  sharpness  of  the  rise  of  this  gust  is  governed  by 
-m  and  the  decrease  by  ~n,  since  its  curve  approaches 
that  of  ( \~e~ml )  near  the  origin  and  finally  becomes 
asymptotic  to  e~nl. 

In  the  case  of  a  rotating  gust  it  is  probably  more 
logical  to  use  the  transient  forms  that  represent  the 
gust  as  disappearing  in  time  instead  of  being  persistent. 

CONTROLLED  MOTIONS 

When  considering  controlled  motions,  it  is  often  just 
as  reasonable  to  assume  that  the  airplane  is  under  a 


THE  CALCULATION  OF  DISTURBED  MOTIONS  OF  AN  AIRPLANE 


kinematic  constraint,  or  prescribed  acceleration,  im¬ 
posed  by  the  control  as  to  assume  that  the  pilot  uses 
the  control  in  an  arbitrary  way.  This  assumption 
leads  to  the  inversion  of  the  integration  problems  here¬ 
tofore  considered,  because  the  motion  of  the  airplane 
is  itself  predetermined  and  the  forces  and  moments  (or, 
more  properly,  accelerations)  required  to  be  supplied 
by  the  controls  are  calculated  by  differentiation.  The 
ability  of  various  control  devices  to  produce  a  given 
maneuver  of  the  airplane  may  thus  be  compared  and 
the  degree  of  coordination  required  of  the  other  con¬ 
trols  may  be  studied. 

The  foregoing  procedure  is  a  particularly  useful  way 
of  studying  the  lateral-control  effectiveness  in  turns. 
Turn  maneuvers,  which  usually  begin  and  end  in  level 


t 

(C)  v0=  Ke-n'(l-e~mi) 


Figure  4.— Curves  of  different  formulas  for  representing  gusts. 

flight,  may  be  described  by  means  of  a  few  sine  or 
cosine  terms.  For  example,  the  angle  of  bank  #  may 
be  given  by 

tp= constant -\-Ai  cos  nt-\-A2  cos  2nt-\-etc.  (27) 
(See  fig.  5.)  The  rate  of  rolling  at  every  instant 

63689-  '36 - 2 


( 

naturally  follows  by  differentiating  this  equation.  If 
the  turn  is  to  be  “perfect,”  that  is,  with  no  sideslipping, 
the  rate  of  yawing  throughout  must  bear  a  definite 
relation  to  the  angle  of  bank,  namely, 

r—yj  sin  ip  (28) 

L'o 

or,  simply, 

r=Yr<p  if  (p  is  under  30° 

Differentiating  the  expressions  for  p  and  r  gives  the 
accelerations  in  rolling  and  yawing  and  hence  the 


Figure  5. — Specifications  for  a  turn  maneuver  in  which  the  constraints  are  given  by 
v-VAAi+Ai  cos  nt+  V\A\  cos  2 vt;  j f=jj^  <p. 


moments,  which  will  arise  from  two  sources:  the 
reactions  due  to  natural  stability  and  the  reactions 
produced  by  the  displaced  controls.  The  reactions 
arising  from  the  motions  are  found  by  combining  the 
known  stability  derivatives  with  the  angular  velocities 
p  and  r,  obtained  from  the  specification  equations  (27) 
and  (28).  The  parts  of  the  moments  necessarily 
supplied  by  the  controls  are  then  obtained  by  deduct¬ 
ing  these  from  the  total  moments.  In  the  case  of  the 
aileron  control,  secondary  moments  in  yaw  result  from 
the  application  of  rolling  moment,  which  modify  the 
amount  of  rudder  control  displacement  necessary. 


CONTROL  AGAINST  GUSTS  OR  ENGINE  FAILURE 

In  order  to  deal  with  attempted  control  of  a  given 
disturbance  it  is  important  to  consider  that  there  is 
invariably  a  lag  in  the  pilot’s  reaction  in  countering 
the  motion.  In  these  cases  it  is  possible  to  assume 
that  the  disturbance  arises  instantly,  or  nearly  so 
(whether  persistent  or  not),  and  that  the  pilot’s  dis¬ 
placement  of  the  corrective  control  takes  place  accord¬ 
ing  to  the  law 

5— 50(1  —e~nt)  (29) 

(see  fig.  4(B))  where  50  is  the  assumed  maximum  con¬ 
trol  deflection,  which  occurs  more  or  less  quickly  as 
— n  is  made  large  or  small. 
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EXAMPLES  SHOWING  APPLICATION  OF  OPERATIONAL 
METHODS  TO  PROBLEMS  OF  AIRPLANE  MOTION 

The  following  examples  illustrate  the  application  of 
the  various  methods  to  specific  problems  of  airplane 
motion.  The  airplane  assumed  in  these  calculations 
is  a  typical  2-passenger  machine  having  the  following 
c-1  lar  act-eristics: 

CHARACTERISTICS  OF  TYPICAL  AIRPLANE 

Type:  Monoplane;  aspect  ratio  6;  rectangular,  rounded 
tip,  Clark  Y  wing;  dihedral  angle,  1°. 
Dimensions: 

Gross  weight _  1,600  lb. 

Wing  span _  32  ft. 

Wing  area _  171  sq.  ft. 

m£x2 _  1,216  slug-ft.2 

mkz _  1,700  slug-ft.2 

Stability  derivatives: 


The  terms  L{)  and  Ar0  represent  the  accelerations  due 
to  the  constant  control  moments  suddenly  applied  at 
/— 0.  They  are 

T  C,qSb 
L°~  mtcx 

(31) 

N=Cd §b 

0  mkz 2 

The  substitution  of  D  for  d/dt,  and  the  rearrange¬ 
ment  of  the  equations  result  in 

Dh — (jp  +  D  UQr = 0 

— Lvv+(D  LqJ'P  Lrr=L0  <  (32) 

— Nvv  ~  Npj>  +  (D— Nr)  r = Ar0  , 

Since  the  rolling  motion  is  desired,  the  equations  will 
be  solved  for  p.  The  algebraic  solution  is: 


Uo 

(ft  ./sec.) 

Cl 

Lv 

hr 

UoLt 

Nr 

UqNv 

150 

0. 35 

-5. 44 

1.  11 

-2.16 

-0.207 

-0.  913 

5.52 

88.5 

1.0 

-3.  23 

1.88 

!  -1.11 

-.301 

-.  663 

2.  04 

06 

*1.8 

-2.  46 

2.51 

;  -1.66 

-.310  , 

-.977 

1.  46 

D 2 

0 

DUq 

—L, 

L0 

~Lr 

-N„ 

Vo 

C D-Nr ) 

D 2 

-(1 

DU0 

—L, 

(D—Lp) 

-Lr 

-N, 

-NP 

(D-Nr) 

The  calculated  principal  lateral-stability  derivatives 
of  this  machine  given  with  the  other  characteristics 
refer  to  motions  of  a  set  of  axes  fixed  in  the  airplane 
but  so  inclined  that  the  X  axis  points  in  the  direction 
of  the  relative  wind  in  straight  flight  at  the  lift  coeffi¬ 
cient  specified.  The  axes,  nevertheless,  move  with  the 
machine  during  the  small  oscillations  considered  and 
hence  depart  slightly  from  instantaneous  reference  axes 
fixed  in  the  wind  direction. 

ILLUSTRATION  OF  SOLUTION  WITH  CONSTANT  DISTURBANCE 

TERM 

Example  I,  Rolling  motion  produced  by  deflecting 
ailerons  at  low  speed: 

(a)  Assume  the  machine  to  be  in  level  steady  flight 
at  a  speed  of  88.5  feet  per  second  (CL— 1.0)  and  that 
a  rolling  moment  corresponding  to  Ci  — 0.04,  with  an 
adverse  yawing  moment  Cn=  —  0.01,  is  applied  sud¬ 
denly  at  the  time  t  —  0.  This  condition  corresponds 
approximately  to  a  full  deflection  of  ordinary  ailerons 
at  this  speed. 

(b)  The  equations  of  the  motion  in  the  three  degrees 
of  lateral  freedom  may  be  set  up  without  including  the 
expressions  for  the  lateral  air  force,  since  this  force  is 
small  and  may  be  neglected  in  this  case.  The  equa¬ 
tions  arc: 

dv  TT 

dT^-rU* 

-^—vLt-\-'pLp-\-rLT-\-La  ■  (30) 


which  is  then  reduced  to  the  form  required  for  expan¬ 
sion  in  equation  (6), 

,_M)r  jmN  (34) 

P~F(D)L°+F(D)JS° 

The  calculation  of  the  various  polynomials  in  D  results 
in: 

h  (D)  =  TF-NJ.P+  U0NVD 
fi  (D) =LrD~—  U„L  VD 

F(D)  =N-  (L,+Nr)D>+  (lPNr-lrNp+  U0NV)D*‘  (35) 
+  u(LvNv-LvN--fjL))T)+g{LvN-LrNt) 


<j-=vNv+VNv+rNr+N0 


At  the  assumed  speed  of  88.5  feet  per  second 
(CY=  1  -0)?  the  constant  rolling  and  yawing  accelera¬ 
tions  are 

1.68 

(36) 

A0=— 0.301 

Using  these  numerical  values,  combined  with  those 
given  for  the  stability  derivatives,  the  polynomials  in 
D  become 

ji  (D)  —  Dz  +  0 . 663D2 + 2 . 04D 

f2(D)  =  1.8SD2+l.UD  (35a) 

F(D)  =D4+3.89D3+4.75D2+10.33D— 1.13 

In  order  to  perform  the  expansion  of  by  Heavi¬ 

side’s  theorem  it  is  necessary  to  determine  the  roots  of 
the  complementary  equation  F(D)  =  0.  When  the 
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polynomial  F(D)  is  plotted  as  a  function  of  a  real  vari¬ 
able  (Z>),  two  real  roots  of  this  equation  are  found: 


Xj=  — 3.41  | 
X2=  0.104  I 


(37) 


Bv  the  use  of  vector  diagrams  (see  fig.  2)  and  the  plot¬ 
ting  of  a  map  of  the  polynomial  considered  as  a  func¬ 
tion  of  a  complex  variable  (D=x  +  iy),  the  following 
root  was  found  by  interpolation: 

X3  — 1.78  (cos  1.73  +i  sin  1.73)  (38) 

An  additional  complex  root  that  is  the  conjugate  of  X3 
is  known  to  exist  and  completes  the  four  roots  of  the 
fourth-degree  equation, 

X4=L78  (cos  1.73 -i  sin  1.73)  (39) 

The  next  step  is  to  set  up  the  integration  equation 
and  perform  the  indicated  operations.  Since  the 
applied  control  moments  L0  and  AT0  are  constants,  form 
(6)  will  be  used 


The  logarithm  of  the  result  naturally  occurs  in  the  form 
x+iy.  Plotting  this  point  on  the  paper  and  construct¬ 
ing  from  it  a  line  parallel  to  X3,  we  obtain  the  locus  of 

ino./i  ( x3)  jo+/2  ( X3) 

for  various  values  of  t  (see  fig.  3).  The  angles  of 
the  final  points  are  given  by  the  ordinates  of  these 
logarithms  and  the  absolute  lengths  by  the  antiloga- 
rit-hms  of  the  abscissas.  The  final  points  are  found  to 
lie  on  a  logarithmic  spiral  whose  radius  decreases  with 
the  time  (time  measured  as  angle)  showing  the  damping 
of  this  component  of  the  motion.  The  summation 
over  the  two.  conjugate  roots  X3  and  X4  is  accomplished 
without  any  further  calculation  by  merely  doubling 
the  abscissas  of  the  points  plotted  above,  as  has  been 
pointed  out.  The  values  thus  obtained  are  listed  in 
the  following  table: 


n  MD)r  I  MD)at  r 

y~~F(D)L^ iw  0  01 


m 

my 


.JjM 


+N, 


M o) 

F(0)' 


■E 


AW_e 

\F'(\y 


(40) 


The  various  terms  to  be  substituted  in  this  formula 
are  found  to  be: 


fi  (o)  j  JM 

F(0r°^F(0) 


A70=0 


/1(X)io+/2(X)A^o=1.68X3+0.54X2+3.09X 


(41) 


XA'(X)=4X4+ 11.67  X3+ 9. 49X2+ 10.33  Xj 
These  terms  are  to  be  calculated  for  the  four  (real  and 
complex)  values  of  the  roots.  In  the  case  of  the  real 
roots  the  calculation  is  made  without  resorting  to 
graphical  methods.  For  Xx=— 3.41,  the  value 

^(Xl)£++l)iVo=  -0.484  results, 

Xi^T  (Xi) 

and  for  X2—0.104 

/l(X2)£o+/2(X2)iVo_n  07>7 

X2F'(X2)  0,277 

It  will  be  convenient  to  perform  graphical  calcula¬ 
tions  to  determine  the  other  parts  of  the  solution, 
corresponding  to  the  complex  terms.  This  result  is 
accomplished  by  calculating  the  square,  cube,  and 
fourth  power  of  the  absolute  length  of  X3  and  by  mul¬ 
tiplying  each  of  these  values  by  the  proper  coefficients 
in  the  polynomials  j{D)  and  F(D ).  By  vector  addition 
the  value  of  the  first  polynomial  was  determined  as 
jx ( X3) L0 +/2 (X3)  2V0 — 3 . 99  (cos  5.23+i  sin  5.23)  (44) 


(42) 

(43) 


Table  of  values  obtained  from  graphical  construction 


t 

(Seconds) 

y+,(X)Lo+MA)Wo,xt 

A  AF'(A) 

(For  X3  and  X«) 

0 

0. 184 

.2 

.186 

.4 

.  166 

.6 

.  130 

.8 

.082 

1.0 

-.030 

1.5 

-.080 

2.0 

-.110 

At  the  time  t— 0,  eM  will  be  unity  so  that  the  initial 
condition  of  zero  rate  of  rolling  should  be  given  by  the 
sum  of  its  coefficients.  The  summation 

—  0.484  +  0.277  +  0.184  —  —0.023 

shows  how  nearly  this  condition  is  attained.  Figure  6 
shows  the  resultant  rate  of  rolling  and  the  components 
of  the  solution  corresponding  to  each  of  the  four  roots, 
X*.  In  addition  to  the  rolling  curve  obtained  by  the 
foregoing  methods,  other  curves  obtained  by  step-by- 
step  integrations  of  the  same  equations  of  motion  are 
given.  In  the  calculation  of  these  curves,  steps  of  one- 
tenth  and  one-twentieth  second  were  taken,  which 
resulted  in  the  differences  shown. 


and  the  second 

X3F/(X3)  =  40.6  (cos  5.60  +  i  sin  5.60).  (45) 

Since  the  quotient  of  these  values  is  to  be  multiplied 
into  <+*  for  a  series  of  values  of  t  it  will  be  convenient 
to  use  the  logarithm  of  this  quotient,  simply  adding  to 
it  the  various  values  of  X3£  for  which  the  calculation  is 
to  be  made.  This  logarithm  is 

n  f 1  (  X3)  Lq  +/ 2  (  X3)  Nq 

log  x3F'(x3) 


ILLUSTRATION  OF  SOLUTION  WITH  VARIABLE  DISTURBANCE 
TERMS 

Example  II,  Sideslipping  during  2-control  turn 
maneuver: 

(a)  Assume  the  airplane  to  perform  a  specified  bank¬ 
ing  maneuver  by  application  of  a  variable  rolling  mo¬ 
ment.  If  no  yawing  moments  (from  either  rudder  or 


=  (log  3.99—log  40.6)  +i (5. 23  — 5.60)  —  —2.32  —  0.38  i 


(46) 
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ailerons)  are  applied,  the  natural  stability  of  the  air¬ 
plane  will  cause  it  to  turn  in  a  direction  appropriate  to 
the  direction  banked.  Such  a  turn  is  called  a  “2-con- 
trol  turn/’  inasmuch  as  only  two  (ailerons  and  elevator) 
of  the  three  available  controls  are  used.  Since  there 
will  not  be  a  very  perfect  coordination  between  the 
banking  and  yawing,  some  sideslip  will  result.  It  is  of 
interest  to  know  the  approximate  amount  of  this 


0  t  2 


Time,  seconds 

Figure  6.— Result  of  sample  computation  compared  with  step-hy-sfep  integrations; 
example  I.  Rolling  motion  following  sudden  deflection  of  ailerons.  Typical  1,600. 
pound  airplane.  Cl  =  1.0;  Ct=0.04;  Cn-~ 0.01. 

sideslip  during  such  a  turn  in  studying  the  practica¬ 
bility  of  2-control  operation. 

(b)  The  first  step  in  this  problem  is  the  determination 
of  a  suitable  expression  for  the  banking  part  of  the 
maneuver.  It  was  considered  that  the  pilot  would 
naturally  conform  his  use  of  the  control  to  the  desired 
motion  of  the  airplane  rather  than  move  the  control  in 
a  predetermined  way  and  accept  whatever  motion  of 
the  machine  followed.  Hence  it  seems  more  logical 
to  specify  the  banking  motion  itself  rather  than  to  try 
to  predetermine  a  law  of  application  of  rolling  moment. 

The  airplane  is  thus  assumed  to  be  constrained  in 
banking  by  the  aileron  control  so  as  to  follow  a  well- 
executed  bank  maneuver  and  recovery.  The  usual 
procedure  in  making  a  turn  is  to  bank  the  machine  up 
to  a  definite  angle,  holding  this  angle  steadily  for  a  short 
time  while  in  the  steady  part  of  the  turn,  and  then  to 
recover  to  level  flight  on  the  completion  of  the  desired 
angle  of  turn.  A  curve  representing  such  a  relation  of 
bank  angle  against  time  may  be  represented  by  a  series 
of  only  two  cosine  terms  with  a  constant  defining  the 


initial  and  terminal  conditions  of  level  flight,  or  zero 
bank  angle.  (See  fig.  7.)  For  a  fairly  sharp  turn  with 
this  small  airplane  the  time  required  will  be  about  6 
seconds  if  the  maximum  angle  of  bank  is  30°.  The 
specification  decided  on  is: 

Bank  angle,  <p=0. 327  — 0.262  [cos  cos  2 1]  (47) 

which  reaches  a  steady  value  of  30°,  and  gives  level 
flight  at  t~  0  and  2x  seconds.  The  rate  of  rolling 
is  the  rate  of  change  of  this  angle  of  bank;  or 

p=<7<=0'262  sin  *+0-131  sin  2t  (48) 

A  constraint  of  the  machine  in  one  of  its  degrees  of 
freedom  having  thus  been  specified,  it  is  only  necessary 
to  consider  the  equations  for  free  motion  in  the  remain¬ 
ing  two  degrees.  As  before,  the  lateral  motion  will  be 
assumed  to  be  independent  of  the  longitudinal.  There 
remain  only  the  sideslipping  and  yawing  motions  to 
be  considered.  Their  equations  are: 


dv  TT 

Tt=OV-rUo 

jt=vNv+VNv+rNr 


(49) 


Although  the  equations  contain  the  rate  of  rolling  and 
the  angle  of  bank,  these  are  to  be  considered  as  known 


Figure  7.— Result  of  computation;  example  II.  Sideslip  during  a  2-control  turn 
maneuver. 


from  equations  (29)  and  (30)  and  are,  in  fact,  to  be 
used  as  the  disturbance  terms.  Calling 


y _ _  0 

U-IT« 

N0=pNj, 


(50) 


and  rearranging  the  equations  as  in  the  other  problems: 


v  fj  \ 

Co 

-UttN?-+(D-N,)r=Ka 

L  o 


(51) 
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Solving  algebraically  for  v/U0: 


or 


where 


Y0 

1 

IN 

(. D-N, ) 

D 

! 

1 

-  U0NV 

C D-Nr ) 

v 

Uo 


MD)y  ,  MPj 

~F(D)Xa^F(D) 


No 


f\(D)=D—NT 

MD)=- 1 

F(D) =FP-NrD+  U0NV 


(52) 


(52a) 


(53) 


If  the  airplane  is  to  maintain  its  altitude  while 
turning,  the  speed  must  be  adjusted  to  give  a  higher 
lift  than  that  at  an  equal  lift  coefficient  in  level  flight. 
At  an  assumed  lift  coefficient  of  1  the  speed  necessary 
to  maintain  altitude  while  turning  at  30°  bank  is 
found  to  be  95  feet  per  second.  Actually,  if  this 
speed  is  held  throughout  the  specified  maneuver,  the 
longitudinal  path  will  be  accelerated  somewhat;  this 
condition  will  be  neglected  in  the  present  problem. 
The  necessary  stability  derivatives  calculated  for  the 
new  condition  are: 

1V,=  -0.712 

UoN„=2AO  (54) 

N,=  -  0.323  | 

The  “disturbance  effects”  Fn  and  Nn  are  (see  equations 
(46),  (47),  and  (49)) 


7o=0. 111-0.0888  cos  < — 0.0222  cos  2 1 

(55) 

N0= -0.0846  sin  t  -0.0423  sin  2 1 
and,  finally, 

c°s  i 

—  0.0222^|^j  cos  2t— 0.0846p|^j  sin  t  (56) 

—  0.0423  sin  2 1 

For  the  expansion  of  these  terms  in  the  integration 
equations  (6),  (8),  and  (9),  it  is  necessary  to  know  the 
roots  of  F{D)  —  0.  These  are 

%  Y±W2-4£/„fV„ 

2  (57) 

X=-0.35G±1.51i  ) 

Since  both  these  roots  are  complex,  the  operations 
indicated  in  the  integration  equations  were  performed 
graphically  in  the  manner  previously  shown. 

The  results  of  these  calculations  are  shown  in  figure  7 . 
The  fact  that  the  error  in  meeting  the  zero  sideslip 
condition  at  the  start  of  the  maneuver  was  very  small 
(even  though  the  graphical  construction  of  several 
terms  was  required)  gives  an  indication  of  the  accuracy 
of  the  calculation. 


Langley  Memorial  Aeronautical  Laboratory, 
National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  February  19 ,  1936. 


APPENDIX 


EVALUATION  OF  ELEMENTARY  OPERATORS 

A  simple  differential  equation  may  be  used  to  illus¬ 
trate  briefly  Heaviside’s  method  of  evaluating  more 
elementary  operational  forms.  Consider  tlie  case  of  an 
airplane  executing  pure  rolling  motion  under  the 
influence  of  a  suddenly  applied  rolling  moment  of 
magnitude  mkx2La>  which  produces  the  impulsive  accel¬ 
eration  La  in  roll.  The  equation  of  motion  may  be 
written: 

| \=PLP+L0  (58) 


in  which  both  p  and  L„  are  supposed  to  have  the  value 
zero  at  the  time  £=0. 

The  solution  of  this  equation  as  ordinarily  found  will 
consist  of  two  parts,  one  of  which  is  a  solution  of 

1-^=0  <59> 

the  “complementary  equation. 77  In  effect,  Heaviside 
wrote  both  equations,  (58)  and  (59),  as  one  by  intro¬ 
ducing  a  discontinuous  function  of  t  into  (58).  Thus, 
(substituting  the  usual  D) 

Dp-Lpp=l(t)L0  (60) 

where  the  symbol  1(0  is  termed  the  “unit  function/7 
and  is  supposed  to  have  the  value  zero  until  the  time 
t= 0  and  to  take  the  value  1  thereafter.  The  algebraic 
solution  of  (60)  is  then  written 


(D—Lp)~l  1(f)  —  (eLv' —  1 )  (65) 

J-jp 

The  final  solution  of  the  original  equation  (1)  follows  as 

p=l(<)f°(eV-l)  (66) 

-LJp 

Such  forms  as  the  left  side  of  equation  (60),  involving 
the  symbol  D ,  are  termed  “operators.77  Equations 
(6)  to  (9)  of  the  text  are  to  be  considered  as  evaluations 
of  the  more  complex  operators  f(D)/F(D)  along  the 
above-indicated  lines.  The  evaluation  of  a  number 
of  such  forms  is  given  in  reference  5. 

Equation  (6)  of  the  text  is  a  shorthand  method  of 
arriving  at  the  foregoing  solution.  For  the  present 
problem  this  formula  is: 

—l(t)  ^~L  (67) 

p”F(D)1WjL°“1W  F{b)L^^\F'{\)  1  ; 


and  the  various  terms  are: 


f(D)  =  l 
F(D)=D—LP 
jf(0)  =  l 
F{0)  =  -L, 


and  it  is  required  to  evaluate  the  form 

zJ=rP1(0 

The  procedure  is  to  expand  the  fraction  by  the  bi¬ 
nomial  theorem  in  ascending  powers  of  L„,  thus, 

(D-Lp)-'=D-'+D-2LP+D-3LpS+  .  +etc.  (61) 

Since 

f  (62) 

D-21{t)==ffiit)dtdt=i(t)-;  etc., 
performing  the  indicated  integrations  results  in 

(D-Lp)-'l  (0  =  1  (<)(*+; |t+^T+  •  •  •  +  etc')  (63) 

If  this  series  is  multiplied  throughout  by  Lp  it  becomes 
identically  the  series  for  eL ^  except  for  the  term  1, 
that  is 

[  LpiD-Lp)-1  + 1]1  (0  =  1  (t)e**  (64) 


m= 1 

The  substitution  of  these  terms  in  (67)  results  in 


as  before. 
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Positive  directions  of  axes  and  angles  (forces  and  moments)  are  shown  by  arrows 


Axis 

Force 
(parallel 
to  axis) 
symbol 

Moment  about  axis 

Angle 

Velocities 

Designation 

Sym¬ 

bol 

Designation 

Positive 

direction 

Designa-  Sym- 
tion  bol 

Linear 
(compo¬ 
nent  along 
axis) 

Angular 

Longitudinal _ 

X 

X 

Rolling—  L 

Y - >Z 

Roll _  0 

u 

V 

Lateral.  _  __  __ 

Y 

i  Y 

Pitching -  M 

Z - >X 

Pitch _  6 

V 

q 

Normal _ 

Z 

Z 

Yawing—  N 

X - >Y 

Yaw _ 

IV 

r 

Absolute  coefficients  of  moment  Angle  of  set  of  control  surface  (relative  to  neutral 

i  n  position),  5.  (Indicate  surface  by  proper  subscript.) 

Cl~qbS  Cm~qcS  Cn~([bS 

(rolling)  (pitching)  (yawing) 

4.  PROPELLER  SYMBOLS 


D,  Diameter 

pf  Geometric  pitch 

pjD,  Pitch  ratio 

Vr,  Inflow  velocity 

V9J  Slipstream  velocity 

T 

T,  Thrust,  absolute  coefficient  =  j)* 

Q,  Torque,  absolute  coefficient  CQ  = 


P, 

Cs> 

y, 
n , 


Power,  absolute  coefficient  CP  = 

Speed-power  coefficient  = 
Efficiency 

Revolutions  per  second,  r.p.s. 


Effective  helix  angle  =  tan  1 


5.  NUMERICAL  RELATIONS 

1  hp.  =  76.04  kg-m/s  =  550  ft-lb./sec.  1  lb.  =  0.4536  kg. 

1  metric  horsepower  =  1.0132  hp.  1  kg  =  2.2046  lb. 

1  m.p.h.  =  0.4470  m.p.s,  1  mi.  =  1,609.35  m  =  5,280  ft. 

1  m.p.s.  =  2.2369  m.p.h.  1  in  =  3.2808  ft. 


